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ABSTRACT 

We determine the contribution of AGN to the mid-IR emission of luminous infrared galaxies (LIRGs) 
ai z > 0.6 by measuring the mid-IR dust continuum slope of 20,039 mid-IR sources. The 24 sources 
are selected from a Spitzer/MIPS survey of the NOAO Deep Wide-Field Survey Bootes field and have 
corresponding 8 /xm data from the IRAC Shallow Survey. There is a clear bimodal distribution in 
the 24 /im to 8 /xm flux ratio. The X-ray detected sources fall within the peak corresponding to 
a flat spectrum in vi^^, implying that it is populated by AGN-dominated LIRGs, whereas the peak 
corresponding to a higher 24 /im to 8 /xm flux ratio is likely due to LIRGs whose infrared emission 
is powered by starbursts. The 24 /im emission is increasingly dominated by AGN at higher 24 /im 
flux densities (f24): the AGN fraction of the z > 0.6 sources increases from 9% at f24 ~0.35 mJy to 
74±20% at f24 « 3 mJy in good agreement with model predictions. Deep 24 /im, small area surveys, 
like GOODS, will be strongly dominated by starburst galaxies. AGN are responsible for ~3-7% of 
the total 24 /im background. 

Subject headings: galaxies: active — galaxies: starburst — X-rays — infrared: galaxies 
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1. INTRODUCTION 

Luminous and Ultraluminous infrared galax- 
ies (LIRGs) have huge infrared luminosities 
(LiR > 10 Lq), corresponding to an extremely 
active phase of dust enshrouded star formation and/or 
AGN activity. They are under- luminous at rest-frame 
UV wavelengths because they are reprocessing and 
radiating much of their energy in the infrared (e.g., 
iSanders k Mirabell [199^ . LIRGs become an mcreas- 
ingly significant pop ulation at high redshifts (e.g., 
iLe Floc'h et al.l I2004D and are likely to dominate the 
luminosity density at z > 1, representing an important 
phase in the build-up of massive galaxy bulges and in 
the growth of their central super-massive black holes 
(SMBHs). To understand these processes, it is essential 
to separate the contribution of AGN and starbursts to 
the IR luminosity of LIRGs. 

Recent work has shown how the IRAC color-color di- 
agram and MIPS 24 /im to 8 /xm c olor can be used 
to identify AGN-dominated LIRGs iLacv et all 120041: 
ISaiinaet a,ljl2005t IStern et a,l ]|20M I Ya,n et a,lJl2004V In 
this paper we demonstrate how the 24 /xm to 8 /im flux 
ratio {C = logio[i^/,y(24/im)/z^/jy(8/im)]) can be used to 
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disentangle the contribution of AGN and starbursts to 
the total reprocessed mid-IR (« 5 — 25 /im) emission as 
a function of 24 /im flux. The utility of the mid-IR slope 
as an AGN/star burst diagnostic wa s first demonstrated 
on samples of nearby ffi^ 5' galaxies ijSanders et alJlT98^ 
Soifer fc Neugebaueil I199H) and further developed by 
Laurent et alJ l|2000() . The powerful UV and X-ray emis- 



sion from an AGN can heat the surrounding dust to a 
range of temperatures up the sublimation temperature of 
«1500 K, resulting in a shallow mid-IR power-law contin- 
mim re.g..lErvis et a,lJl 994lKlaas et a]J200lULutz et alJ 
119981) . In contrast, hot stellar sources in star- forming 
galaxies do not typically heat dust to these temperatures, 
resulting in a steeper mid-IR continuum slope. If the ob- 
served 24 /im and 8 /im flux densities are not heavily 
influenced by absorption and emission features, AGN- 
dominated LIRGs should have smaller values of C than 
starburst-dominated LIRGs. Cooler dust temperatures 
which are characterized by steeper mid-IR continua may 
also be expected in sources containing only a weak or 
heavily obscured AGN. We will return to these issues in 
Section^ In this work, we use the mid-IR slope as an 
AGN/star burst diagnostic applied to a large statistical 
sample of sources observed with the Spitzer Space Tele- 
scope. 

2. DATA 

The sample for this study was drawn from Spitzer ob- 
servations of th e NOAO Deep Wide-F ield Survey Bootes 
field fNDWFS: IJannuzi & DevllT999l) . The 9.93 deg^ re- 
gion has been mapped by Spitzer/MIPS at 24 /xm (f24; 
comprising 28,391 sources down to a 5a depth of ^0.3 
mJy; version made on July 8th 2005). An 8.5 deg^ region 
of the field has been mapped in all four S pitzer /IRAC 
band s by the IRAC Shallow Survey (Eisen liardt et all 
|2004() . reaching 5a depths at 8 /im of «0.076 mJy (for 
the fluxes measured in 6" apertures). The 24 fim cata- 
log was created by running a PSF fitting technique us- 
ing DAOPHOT. In this paper, only the 25,681 24 /im 
sources whose positions overlapped with data in both 
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Fig. 1. — Histograms showing (a) f for all 24 /^m sources, and (b) all 24 /im sources excluding those with z < 0.6, split into £24 bins. 
The histograms are for £24 > 0.3 mjy (all sources; solid black line), 0.3 < £24 < 0.5 mjy (short-dashed blue line), 0.5 < £24 < 0.7 mjy 
(dot-dashed green line), 0.7 < £24 < 1.2 mjy (dot-dot-dot-dashed red line), and £24 > 1.2 mjy (long-dashed purple line). Also plotted is the 
f distribution o£ X-ray sources detected in the XBootes survey (dotted orange line). The number o£ sources contributing to each histogram 
is listed in the figure. The shaded histograms show the distribution o£ the 24 fira sources with no 8 fim counterparts assuming an 8 fim 
flux density limit o£ 3 times the locally calculated sky RMS. 



the NDWFS and IRAC Shallow Survey are considered 
(w 8.2 deg^). Optical and infrared identifications were 
determined by matching the 24 /im catalog to the multi- 
wavelength catalogs using a 2" matching radius. We used 
the 8 /im flux densities (fg) measured in 6" apertures cor- 
rected to the total flux density assuming a point source 
profile. In cases with no 8 /Ltm catalog detection (1588 
sources), fg was measured directly from the images; when 
no signal was found (557 sources), the fg limit (defined 
as 3 times the locally measured sky RMS), was used. In 
practice, the inclusion or exclusion of these sources does 
not significantly affect any of our conclusions. 

The Bootes region has also been surve yed by the 5- 
ks Chandra X- ray Observatory ( XBootes: iMurrav et all 



I2m5l iKenter et al1l2(i05l: iBrand et al.' 2006). The AGN 
and G alaxy Evolution Sur vey (AGES; Kochanek et al. in 
prep.; iBrown et al.ll2005|) has also targeted the Bootes 
field, providing spectroscopic redshifts for 98% of the 
hi > 0.3 mJy 24 /zm sources with / < 21.5 (optical 
point sources) and f24 > 0.5 mJy sources with / < 20.0 
(optical extended sources). As described in Section 3, we 
use this survey to exclude z <0.6 sources in which strong 
PAH emission may fall into the observed 8 /im bandpass. 

3. ANALYSIS 

We calculated C for each 24 /im source as a crude mea- 
sure of the mid-IR continuum slope. This provides an 
estimate of the relative luminosities of the hot (« 750 
K at z wl) and warm (« 250 K at z «1) dust compo- 
nents. Histograms of ( for different f24 bins are plotted 
in Figure^. The distribution has 4 peaks, located at 
( » -1.5, C ~ -0.4, C ~ 0, and C ~ 0.5. As we discuss 
below, we identify these peaks with stars, low-redshift 
normal star-forming galaxies, AGN, and higher redshift 
starburst galaxies respectively. In addition, sources have 
progressively lower values of C at higher f24. Although 
the sources with no 8 /im detection become a larger frac- 
tion at smaller f24, they should only affect the tail of the 
distributions at C > 1-0 and will not significantly affect 
the location or width of either peak. 

3.1. Low redshift contaminants 



Strong PAH emission and silicate absorption can bias 
the broad-band flux density measurements of LIRGs. 
In low redshift (z < 0.3) galaxies, strong 7.7 /im PAH 
emission can bias ^ to smaller values, rendering normal 
star-for ming galaxies and AGN difficult to separate (e.g., 
lYan et aL.,2Q04) . To avoid this bias, we use existing spec- 
troscopic redshifts from AGES and conservatively ex- 
clude the 5642 sources with z < 0.6 from our sample. 
Although AGES is w98% complete for optically bright 
sources, optically faint (/ > 20.0) low-redshift (z < 0.6) 
galaxies may remain in our sample. Figure shows 
that after exclusion of the AGES z < 0.6 sources, the 
C w —1.5 and C ~ —0.4 peaks become only a small resid- 
ual and the distribution of the remaining 20,039 sources 
becomes largely bimodal. By fitting the ^ peaks with 
Gaussian profiles, we estimate the residual contamina- 
tion by low-redshift galaxies in our z > 0.6 sample to be 
only «2%. 

3.2. The C values of X-ray loud AGN 

X-ray data show that AGN-dominated sources pop- 
ulate the C ~ peak. Of the 20,039 24 /im sources 
considered in this work, most of which are likely at 
z > 0.6, 1095 («5%) have an X-ray counterpart in the 
Chandra XBootes survey; their C distribution is over- 
plotted in Figure ^ These sources occupy a narrow 
distribution around C « 0, i.e., a flat spectrum in i^f^ 
consistent with the typical slope of unobscured AGN 
(e.g., Elvis ct al. 1994). Because the XBootes survey 
has a relatively shallow fl ux density limit (fo.5_7 keV ~ 
7.8xl0~i5 ergs cm'^ s~^: lKenter"er al. 2005.) . all 24 um 
sources at z > 0.6 with an X-ray counterpart must con- 
tain a powerful (L^ > lO'^^ergs s~^) AGN. This does not 
mean that the mid-IR emission from 24 /im sources with- 
out X-ray counterparts cannot be dominated by an AGN. 
The X-ray emission from the AGN could be significantly 
attenuated in LIRGs, even in the hard X-ray band. In 
addition, the X-ray emission from AGN has been shown 
to be highly variable with time (e.g.,[Paolillo ct al. 2004). 
Inclusion of the X-ray data serves purely to demonstrate 
that the ( distribution of 24 /im sources hosting a pow- 
erful AGN peaks at C «0. 
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Fig. 2. — Histograms showing for all 24 sources with 
spectroscopic redshifts z > 0.6 from AGES which show signs of 
AGN activity in their optical spectra (solid red line), which show 
no signs of AGN activity in their optical spectra (dashed blue line), 
and which show both evidence of AGN activity and stellar absorp- 
tion features in their optical spectra (solid green histogram). The 
number of sources contributing to each histogram is listed in the 
figure. 

3.3. The C values of sources with optical spectroscopy 

Of the 20,039 24 /ini sources at z > 0.6, we have optical 
spectra for 1702 opticaUy bright sources from AGES and 
174 opticaUy fainter sources from KECK/DEIMOS. For 
aU optical spectra for which the signal-to-noise ratio was 
sufficient to perform the analysis, we classified the AGES 
sources by eye into those with signs of AGN activity (i.e., 
clear broad line emission lines or high ionization lines), 
only galaxy signatures, or both. Figure |21 shows that the 
sources with AGN signatures fall in a narrow distribution 
around C «0 where we expect AGN-dominated sources. 
The C distribution of sources with no AGN signatures in 
their optical spectra peaks at ( «0.5. There is a smaller 
tail at higher C values than for the general population. 
This is presumably because the AGES sources are all op- 
tically bright and hence do not include the most dusty, 
obscured and/or high redshift starburst galaxies which 
tend to have the highest ( values. The sources with both 
AGN and galaxy signatures have a large range of C val- 
ues, suggesting that they are powered by a combination 
of AGN and starburst activity. 

Figure 131 shows C as a function of redshift for all 174 
z > 0.6 sources for which we have spectroscopic redshifts 
from KECK/DEIMOS and for aU 616 optically bright 
AGES sources which have X-ray counterparts from the 
XBootes survey. Overplotted are LIRG and ULIRG tem- 
plates whose infrared emission is known to be dominated 
by either a starburst (Arp 220, M82, NGC 6090) or an 
AGN (NGC 1068, Mrk 231). The templates show that we 
expect starburst-dominated LIRGs to have significantly 
higher values of C (with ^ « 0.5) than AGN-dominated 
LIRGs at all redshifts except z '--^ 1.1 — 1.6. The pro- 
nounced dip in of starburst-dominated LIRG templates 
is due to a strong 9.7 /im silicate absorption feature enter- 
ing the observed frame 24 /zm band. The possible effects 
of this are discussed in Section ^ The X-ray sources 
span a large redshift range and have C values that are, 
in almost all cases, consistent with the AGN-dominated 
templates. The sources with high ( tend to have lower 
redshifts, perhaps suggesting that the majority of z > 1.2 
starburst-dominated LIRGs are too faint to be included 
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Fig. 3. — = logio[!^/i/(24/im)/i^/^(8/im)] as a function of 
redshift for all z > 0.6 sources for whicli we have spectroscopic 
redshifts from KECK/DEIMOS (174; filled circles) and for all X- 
ray sources from the optically bright AGES survey (616; empty 
stars). Overplotted are LIRG templates which are known to 
be starburs t-dominated (soUd Unes: Arp 220, NGC 6090, M82; 
Silva et alj[T9 98'). and AGN-do minated (dashed h ues: NGC 1068; 
Le Floc'h et a l. 2001, Mrk 2 31; lEgami et al.ll200l averaged Syl/2 
SEDs : IChatzichristouilflfKl) . 



in our 24 fj,m sample (see, e.g.. lGruDDioni et alJl2005 l fig- 
ure 5.), have lower values of ( due to sihcate absorption, 
or have only weak spectral features in the observed-frame 
optical band. 

3.4. The fraction of AGN-dominated LIRGs as a 
function 0//24. 

Assuming that the mid-IR emission from all z > 0.6 
sources in the C ~0 and ^ «0.5 peaks are respectively 
AGN- and starburst-dominated, we can estimate the 
AGN-dominated fraction as a function of f24. The num- 
ber of sources in each peak were estimated by fitting 
the positions, widths, and amplitudes of two Gaussian 
profiles (chosen for simplicity and because they fit the 
data reasonably well). The central position of the AGN 
peak was determined first by fitting a Gaussian profile 
to the C histogram of the X-ray sources. The central 
position of the starburst peak was then determined by 
fitting 2 Gaussian profiles to the C histogram of the 0.3 
mjy < f24 < 0.5 mjy sources while fixing the position 
and width of the AGN peak to their estimated values. 
The widths and amplitudes of the two peaks were then 
simultaneously fit as a function of f24, with their posi- 
tions fixed. Table n lists the resulting fits. The fraction 
of sources whose mid-IR emission is dominated by an 
AGN is estimated from the area under the C ~ Gaus- 
sian divided by the total area. 

Figure^Ji shows that at the brightest f24 , 74±21% of all 
2 > 0.6 sources have their mid-IR emission dominated by 
an AGN. At the faintest f24 probed by our survey (wO.3 
mJy), this fraction decreases to «9%. Thus, in the deep- 
est 24 fim surveys, starbursts will dominate the mid-IR 
emission in the vast majority of sources. We note that 
this is a measurement constrained by the assumption 
that the distributions of the AGN and non-AGN popu- 
lations as a function of ^ can be approximated by Gaus- 
sian distributions. Figure 01d shows the fraction of all 24 
fim sources whose mid-IR emission is dominated by an 
AGN as a function of their 24 /im flux. We assume that 
all z < 0.6 sources are starburst-dominated unless they 
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TABLE 1 

Gaussian fit parameters and AGN fraction for z > 0.6 

SOURCES. 



f24 C = 0.03 peak C = 0.47 peak AGN 

(mjy) Width^ Amp^ Width^ Amp^ fraction 



X-ray 0.17 117.8±4.3 - 



0.3-0.5 


0.24 


98.7±4.5 


0.31 


705.0±5.4 


0.09±0.01 


0.5-0.7 


0.18 


89.6±6.1 


0.25 


251.5±6.0 


0.20±0.03 


0.7-1.2 


0.16 


61.3±4.0 


0.20 


73.3±3.5 


0.38±0.06 


1.2-1.9 


0.15 


31.2±2.9 


0.20 


17.0±1.8 


0.57±0.14 


>1.9 


0.16 


22.9±2.2 


0.32 


4.1±0.9 


0.74±0.21 



the sources with brighter £24 are increasingly hkely to 
populate the C ~ peak. We suggest that this bimodal 
behavior may be due to the presence or absence of a pow- 
erful AGN which heats the surrounding dust to high tem- 
peratures ( ga 1500 K at the dust-sublimation radius; e.g., 
ISanders"et al. 1989,) , resulting in a shallower power-law 
continuum slope for the AGN-dominated sources. Our 
confidence in this interpretation is boosted by the re- 
markably narrow C distribution of the X-ray sources at 

C « 0. 

Making the crude assumption that, on average, the z > 
0.6 sources with higher f24 correspond to sources with 
higher IR luminosities, our results imply that the fraction 
of AGN-dominated sources increases significantly with 
IR luminosity. This is in broad agreement with a num- 
ber of earher studies at low redshift. .Lutz et a l. (199^ 
used ISO to study a sample of brighter, nearby ULIRGs. 
They measured the ratio of the 7.7 /im PAH emission 
feature to the local continuum to determine that while 
only 15% of LIRGs at luminosities below 2xlO^^L0 are 
AGN powered, this number increases to ^^50% at higher 
luminosities. 

Although the ( diagnostic appears to work well in dis- 
tinguishing AGN and starburst dominated sources, we 
note several caveats with its use. Many broad emission 
and absorption line features are kno wn to be present in 
the mid-IR spec trum of LIRGs fe.g.. lHouck et al.ll2005t 
Yan et aDl2005() and this may affect ^ as a function of 
redshift. We have already discussed how the 7.7 /xm PAH 
emission feature depresses C when it passes through the 
8 fim band in z < 0.6 galaxies. At higher z, the most no- 
table broad feature which may affect C, is the silicate ab- 
sorption feature at 9.7/xm. If present, this may strongly 
attenuate the observed 24 /xm emission at zRil.1-1.7, re- 
sulting in a lower C and hence mimicking a shallower con- 
tinuum slope (see Figure |3J. Dusty starburst-dominated 
LIRGs at z R:il.l-1.7 may the refore be mis-classified as 
AGN fe.g.. iDesai etall l2005|) . This should not signifi- 
cantly affect our AGN-fraction determination since, if the 
absorption feature is strong enough, only a small number 
of affected sources will be bright enough to be selected in 
our shallow 24 /im survey. iPearsonI l|2005|) predict «25 
non-AGN sources per square degree with /24 >0.25 mJy 
at z «1. 1-1.7. We therefore expect a maximum of only 
«225 starburst-dominated sources (~1% of our z > 0.6 
sample) to be mis-classified as AGN-dominated sources 
due to the silicate absorption feature. 

It is also possible that an AGN could be heavily em- 
bedded in large amounts of cooler dust and could remain 
undetected in the 8 /xm band even though it dominates 
the bo lom etric infrared emission. lAlonso-Herrcro et alJ 
(2003) and lClavel et all 1)2000 ^ show that nearby Seyfert 
2 galaxies tend to have steeper mid-IR spectra than 
nearby Seyfert 1 galaxies, suggesting that high absorp- 
tion blocks a large fraction of the mid-IR emission from 
the inner torus, at least in low luminosity AGN in the 
local Universe. In the X-ray detected AGN, we notice a 
weak trend of lower C values in the sources with fainter 
optical counterparts, suggesting that obscuration may in- 
fiuence our results. Ideally, one would estimate the con- 
tinuum slope at longer wavelengths to avoid this prob- 
lem. Nevertheless, although these potential effects may 
be important, they do not significantly dilute the bi- 
modal nature of the C distribution and therefore cannot 



^The Icr width of the fitted Gaussian. The uncertainties are <0.01 
"^Amplitude of the fitted Gaussian 

have X-ray counterparts in the XBootes surve y, so this 
calculation should be treate d as a lower limit. IPearsonI 
(|2005j) and iGruppioni et'aLl l)2005|) present phenomeno- 
logical models for the evolution of different populations 
to determine their contributions to the 24 /im source 
counts. Although both models reproduce the observed 
trend of a higher fraction of AGN-dominated LIRGs at 
high f24, the models of .Pearson (.2005.) are in better agree- 
ment with ou r observations. Our re sults agree well with 
the results of iTreister et all (j2005(l who identify AGN- 
dominated sources by their hard X-ray emission for all 
24 /im sources in the GOODS field (see Figure^). 

3.5. The AGN contribution to the 24 /im background. 

To estimate the contribution of AGN to the 24 /im 
background, the z < 0.6 sources were again included, 
and assumed to be starburst-dominated unless they were 
X-ray sources. These results w ere combined wi t h the dif- 
ferential number counts from iPapovich et al.l l)2004|) to 
estimate the contribution of AGN to the 24 /im back- 
ground. Down to the f24=0.3 mJy limit of our sam- 
ple, where we resolve ^ 30% of the 24 /im background, 
AGN contribute «10% of the t otal mid-IR flux. This 
is consistent with the r esults of iFadda et al.l 120021 and 
IFranceschini et al]l2005l To estimate the contribution of 
AGN to the total 24 fim background, our estimates of 
the fraction of sources dominated by AGN were extrapo- 
lated to f24 =0.033 mJy (the flux density hmit of the d ata 
used bv lPapovich et al.ll2004jl . IPapovich et'aLl l|2004|) es- 
timate that they resolve w70% of the background. We 
assume the remaining 30% of the background is com- 
prised of f24 <0.033 mJy starbursts. Assuming that the 
fraction of AGN-dominated sources tends to and 0.1 at 
f 24 =0.033 mJy, we find that AGN contribute 3 and 7% 
of the background respectively. The 24 /xm background 
is dominated by sources with low 24 /im flux (peaking 
at wO.2 mJv: IPapovich et alj[200^ . Although AGN do 
not dominate the mid-IR emission in these sources, they 
could still contribute to a signiflcant fraction of the total 
emission, and therefore make up a larger fraction of the 
total 24 /im background than the above calculation sug- 
gests. In addition, the AGN light could remain obscured 
in the observed 8 /im band but still dominate the emis- 
sion at longer wavelengths. Our result should therefore 
be considered a lower limit. 

4. DISCUSSION 

There is a clear bimodal distribution in the C. values of 
z > 0.6 24 /im sources at C « and C w 0.5. In addition. 
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Fig. 4. — The fraction of (a) all z > 0.6 sources and (b) all 24 /im sources (assuming all z < 0.6 sources are starburst-dominated) whose 
mid-IR emission is dominated by AGN as a function of {24 (filled squares). The vertical error bars are calculated from the uncertainties in 
the Gaussian fit parameters and do not include the uncertainties introduced by the possible inclusion of a small number of optically faint 
z < 0.6 sources or the possibility that some sources may contain a heavily embedded AGN which may be obscured even at 8 fim. The 
horizontal error bars show the uncertaint y on th e mean f24 in each bin. Overplotted are the expected fraction of AGN-d ominated sources 
for the bright-end and burst models of IPearsonI 12003) (dotted and dashed lines respectively) and in (b), the model of 'Gr uppioni et alj 
12305) (long dashed line). Also overplotted in (b) are the results of Treistcr et al. (2005) using a hard X-ray AGN classification in the 
GOODS fields. The points are uncorrected (empty circles) and corrected (empty triangles) for AGN not detected in the X-ray. 



be the dominant effect. 

Obtaining complete samples of optically fainter mid- 
IR sources with spectroscopic redshifts will be important 
in testing the ^ diagnostic and in drawing conclusions 
about the physical nature of high redshift LIRGs. In 
particular, this will allow us to confirm whether the 
fraction of AGN-dominated sources increases with IR 
luminosity at high redshift, and to determine the effects 
of measuring the dust continuum slope at different 
rest-frame wavelengths (can powerful AGN really be so 
embedded that their near-IR emission is absorbed?). 
Redshift information will also be crucial in determining 
the contribution of LIRGs to the star-formation and 
SMBH accretion histories in massive galaxies and in in- 
vestigating whether the AGN- and starburst-dominated 
LIRGS and ULIRGs form an evolutionary sequence. 
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